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Abstract 

The laser cooling of vibrational states of solids has been achieved 
through photoluminescence in rare-earth elements, optical forces in 
optomechanics, and the Brillouin scattering light-sound interaction. 
The net cooling of solids through spontaneous Raman scattering, and 
laser refrigeration of indirect band gap semiconductors, both remain 
unsolved challenges. Here, we analytically show that photonic den¬ 
sity of states (DoS) engineering can address the two fundamental re¬ 
quirements for achieving spontaneous Raman cooling: suppressing the 
dominance of Stokes (heating) transitions, and the enhancement of 
anti-Stokes (cooling) efficiency beyond the natural optical absorption 
of the material. We develop a general model for the DoS modifica¬ 
tion to spontaneous Raman scattering probabilities, and elucidate the 
necessary and minimum condition required for achieving net Raman 
cooling. With a suitably engineered DoS, we establish the enticing 
possibility of refrigeration of intrinsic silicon by annihilating phonons 
from all its Raman-active modes simultaneously, through a single tele¬ 
com wavelength pump. This result points to a highly flexible approach 
for laser cooling of any transparent semiconductor, including indirect 
band gap semiconductors, far away from significant optical absorption, 
band-edge states, excitons, or atomic resonances. 


1 Introduction 

The photon-induced annihilation of thermal quanta from matter, i.e. laser 
cooling, is of key importance in ultracold science and has been instrumental 
in the creation of gas phase quantum condensates, matter-wave interfer¬ 
ometry, and macroscopic tests of quantum mechanics. The laser cooling 
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of solids, in particular, is enabled through photon up-conversion processes 
in which phonons are annihilated from the material m- On this prin¬ 
ciple, bulk cooling of solids has been achieved through photoluminescence 
on specific electronic transitions in Yb-doped solids [1] with demonstrations 
even reaching cryogenic temperatures [SHE]. Photoluminescence cooling, en¬ 
hanced through exciton-phonon coupling, was also employed for the first 
demonstration of laser cooling of an undoped direct band gap semicon¬ 
ductor [9]. However, this method is not effective with indirect band gap 
semiconductors like silicon and germanium due to the need of a phonon for 
momentum conservation and the consequent competition with nonradiative 
recombination. In parallel efforts, efficient narrowband cooling of individual 
phonon modes in solids has also been achieved by means of optical forces in 
geometry-engineered microstructures [Mn] and through Brillouin scatter¬ 
ing m, albeit without net cooling. The search still continues for a universal 
laser refrigeration mechanism that could be applied to any transparent ma¬ 
terial, with any geometry, at any optical wavelength. 

In this context, Raman scattering of light from optical phonons in solids 
has been suggested as a promising alternative technology, especially since 
Raman scattering is available in all materials and also exhibits the requisite 
photon up-conversion i.e. anti-Stokes scattering ismiiiiij. The key distinc¬ 
tion is that Raman cooling involves only virtual excited states while photo¬ 
luminescence mechanisms involve real excitations of the electronic states of 
a system. Photoluminescence thus relies on the existence of suitable natural 
materials, while Raman scattering can be potentially engineered. However, 
two key challenges must be resolved before Raman cooling of a bulk solid can 
be achieved. First, spontaneous Stokes Raman scattering always dominates 
over anti-Stokes scattering in bulk materials due to the event probabilities 
being scaled by uq -|- 1 and no respectively, where no = (^e^o/ksT _ jg 
the Bose-Einstein distribution function. It is essential to invert this imbal¬ 
ance to have any possibility of achieving cooling. Second, the anti-Stokes 
scattering efficiency needs to be enhanced so as to overcome absorption of 
the pump light by the material. This is critical since each anti-Stokes Ra¬ 
man scattering event annihilates fiWQ energy from the material (wo is phonon 
frequency), while each absorption event adds ^photon^o) energy to the 
material. 

Here we show that both of the above requirements can be met in nearly 
any semiconductor by engineering the optical density of states to a suit¬ 
able form. Accordingly, we develop expressions for the spontaneous Raman 
scattering efficiency in a material with a modified photonic density of states 
(DoS), and derive the necessary and minimum condition to achieve cooling 
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by this method. As a specific case, we demonstrate the surprising possibility 
of cooling all the Raman-active phonon modes in undoped crystalline silicon 
using a single telecom wavelength pump. 

Before we proceed, we must recognize the encouraging experimental ef¬ 
forts of past researchers also moving in the direction of this result. It has 
been reported that resonant Raman scattering (using excitonic enhance¬ 
ment) in certain semiconductors [3l |9] can be used to tilt the Stokes vs 
anti-Stokes balance in favor of cooling. Ding and Khurgin showed [3] that 
by pumping GaN at 375.8 nm, anti-Stokes to Stokes ratio of 0.38 can be 
naturally reached. More recently, Zhang et al reported a ratio of 1.43 in 
CdS nanobelts by means of an exciton resonance, firmly establishing that 
anti-Stokes dominance can indeed be achieved. Since this method relies on 
interaction with an exciton, it can only allow narrowband cooling of one or 
a few phonon modes that lie within the resonance, a constraint similar to 
optomechanical cooling [101413] . Additional challenges to this method are 
the large optical absorption that occurs due to strong exciton-photon inter¬ 
action, and its impracticality for cooling indirect band gap semiconductors. 

A different perspective on Raman cooling has been offered by Rand 
using a multi-pump method on specihc electronic transitions of rare earth 
elements like Ce^"*" and Yb^"*". However, as in other efforts each photon 
pumping configuration chosen would target only one phonon mode. 

We propose that engineering the photonic DoS [laiiB] at high-transparency 
wavelengths can be used to invert the imbalance of Stokes vs anti-Stokes 
spontaneous Raman scattering efficiencies, and can also be used to enhance 
anti-Stokes Raman scattering through engineered photonic resonances. This 
proposal is outlined in Fig. This method circumvents the optical absorp¬ 
tion issue confronted by techniques that use excitons mi, can be applied 
to indirect band gap semiconductors, and also avoids band-tail absorption 
confronted in photoluminescence approaches Such DoS engineering 

has been previously employed in experiments with W1 silicon waveguides 
[T8] to enhance spontaneous Stokes Raman scattering by a factor up to 
13. It has also been shown that the naturally occurring photonic DoS of 
a high-Q resonator can be used to enhance anti-Stokes Brillouin scattering 
by orders-of-magnitude while completely suppressing Stokes scattering, thus 
enabling Brillouin cooling m- The possibility of by-design DoS engineering 
to achieve spontaneous anti-Stokes Raman cooling is clear. 
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Figure 1: Concept of density of states (DoS) engineering for Raman cooling of 
an arbitrary number of phonon modes, (a) Stokes scattering (red) dominates anti- 
Stokes scattering (blue) from any chosen phonon mode in bulk media, (b) It 
is proposed that an engineered photonic DoS with a complete band gap can be 
used to (c) suppress Stokes scattering while simultaneously enhancing anti-Stokes 
scattering intensity, as previously demonstrated with Brillouin cooling m- 


2 Raman scattering efficiency with modified DoS 

Consider a spontaneous Raman scattering process with incident photon fre¬ 
quency LJi and scattered photon frequency lo. The scattering probability is 
given by usms] 

Ptt^ 1 

W{uji,uj) = ^uJiUjNi—D{u;)\M\‘^, (1) 

where Ni is the number of incident photons, h is the Planck constant, D{uj) 
is the photonic DoS at the scattered photon frequency, and A4 is the matrix 
element of the scattering process. Here the scattered photon frequency uj 
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can be taken as either the Stokes frequency ujs = coi — ojq or the anti-Stokes 
frequency, loas = where wq is the phonon frequency. If the scattering 

process happens in free space, the DoS is obtained using the linear dispersion 
relation uj = c\k\. 

03 (-) = ( 2 ) 


where k is the wave vector of the photon and c the speed of light. If we 
assume the phonon energy is small compared against the photon energy, we 
can approximate Wj ~ w, and Eqn. Q will recover the well-known W (x uj^ 
relation [16]. On the other hand, in a medium where the photon dispersion 
relation is not given by a; = c|A:|, the scattering probability will be obtained 
by replacing the photonic DoS in Eqn. by the modified DoS. In such a 
generalized case, the scattering efficiency per unit length, 5S, per unit solid 
angle, 5D, is given by 


dS Whuj 

— oc- 

do. Nihui 


oc D{uj). 


( 3 ) 


The major challenge in calculating the Raman scattering efficiency is to 
determine the matrix element Ai in Eqn. 

To quantitatively evaluate the effects of arbitrary photonic DoS on the 
Raman scattering efficiency, we use the formulae found in [201 |2T] for the 
following calculations. The Raman scattering efficiency per unit length and 
solid angle is given by [20] as 


/ Stokes 


/W5\4 Nh 
V C / MoJo 


(l + no)|Rs(^^)P 


for Stokes scattering, and 


=f 

dVl J Anti—Stokes 




\ c J Mojq 




( 4 ) 

( 5 ) 


for anti-Stokes scattering, where N is the number of unit cells per unit 
volume, M is the atomic mass, and no = — 1)“^ is the phonon 

occupation number. R is the Raman tensor element, which depends on the 
scattering angle and the crystal structure of the material. In the simple case 
where the medium is transparent with respect to the pump and scattered 
light, the total Raman scattering efficiencies per unit length, Ss and Sas, 
can be calculated |22| by integrating Eqns. and over the solid angle. 
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Note that the above equations hold only for linear dispersion relations. Fol¬ 
lowing Eqn. if the photonic DoS is modihed Eqns. and become 
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where D(uj,^}) is the modihed DoS at certain solid angle and ^ 3 ( 0 ;) is ob¬ 
tained from Eqn.[^ This DoS ratio factor resembles the Purcell enhancement 
of spontaneous emission m- Note that if the modihed DoS is isotropic, 
D(uj, D) can be replaced by D(a;)/ 47 r simplifying the factor within the inte¬ 
grals to D{uj)/D^{ijj), which is independent of the scattering angle. 

Since the typical phonon energies are small compared against the photon 
energy, the Raman tensor for both processes is symmetrical and has approx¬ 
imately the same values. Erom Eqns. and we see that it is also possible 
to adjust the scattering efficiency in the Stokes and anti-Stokes directions by 
aligning the material’s crystal (Raman selection rules) suitably with respect 
to the anisotropic photonic DoS. 


3 Application to laser cooling 

Now that we have obtained the DoS-modihed Raman scattering efficiencies, 
we can calculate whether this can be used to achieve Raman cooling in solids. 
We consider a highly transparent medium in which Raman scattering occurs 
over a volume of cross-sectional area A and effective length L. We can 
define net scattering efficiencies by rjs = L 8$ and rjAS = LSas-, such that 
Is = Vs Ipump and I AS = VAS Ipump are the scattered intensities, with Ipump 
being the incident light intensity. The experiment should be performed 
at high-transparency wavelengths far away from any electronic resonances, 
exciton resonances, or two photon absorption in the material. R is then fair 
to assume that the only energy exchange mechanisms between the medium 
and the pump laser are residual absorption and Raman scattering. Under 
these assumptions, the net power transferred into the medium is given by 

Pnet — Pabs T Pph,S Pph,ASi (^0) 
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where Pabs is the broad-spectrum absorbed optical power, while Pph,s and 
Pph,AS are the heating and cooling power due to phonon creation and anni¬ 
hilation from Stokes and anti-Stokes processes, respectively. Cooling of the 
entire material will occur when Pnet < 0. 

The absorbed power can be calculated by Pabs = -^Ipumpi^ — 
where a is the absorption coefficient at the incident frequency. Since the 
absorption coefficient is small for a transparent material the absorbed power 
is approximately Pabs = ALa I pump- Presently, we consider only one phonon 
mode for which the scattered optical power in the Stokes sideband is Ps = 
A r]s I pump = ALSs I pump, and similarly Pas = AL Sas I pump for the anti- 
Stokes sideband. For each Raman scattered photon, one phonon is added or 
removed from the system. The heating (cooling) power due to scattering is 
thus obtained by multiplying the creation (annihilation) rate by the phonon 
energy /kuq. 


Pph,S — hwo . — A • L • S sipump, 

hws 


Pph,AS huiQ 


Pas 

huJAS 


= A-L- 


^AS 
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Upon substituting these relations into equation (10), we have the cooling 
condition 


net 


= A-L-1, 


pump 


Wo c, Wo 

a H- Ss - 

W5 WAS 


Sas < 0, 


(13) 


resulting in the necessary and minimum requirement for cooling of the solid 

(14) 


a < 


Wo r, Wo 

- Sas - Ss- 

OJAS ws 


This means that to achieve cooling, the net cooling efficiency per unit length 

exceed the optical absorption coefficient. If the DoS 
is designed such that multiple Raman modes are simultaneously affected, the 
above relation can be suitably expanded for multiple phonon modes into the 
form 


a < 


E 


all modes 


Wo Wo 

- Sas - Ss 

OJAS Ws 


(15) 


4 Raman cooling of silicon 

Amongst common semiconductor materials, intrinsic crystalline silicon is 
extremely transparent in the telecom range |23H25] . and has Raman scat¬ 
tering efficiency comparable to the optical absorption coefficient [2T]. It is 
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thus a promising material with which to demonstrate Raman cooling. In 
this section, we show that the cooling of undoped silicon is practical using 
a telecom pump when a diamond structure 3D photonic crystal is used to 
generate a complete photonic band gap. 

First, we consider a simplification to Eqns. and for cases where an 

isotropic or nearly isotropic photonic DoS is available, which will be used 

later. The ratio of DoS in these equations can then be pulled out from 
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the integral as a simple prefactor D{u!)/{■^^^), where the numerator is the 
photonic DoS of the crystal, and the denominator is the vacuum DoS. We 
then simply need to integrate the variation of the Raman tensor element 
over the solid angle D over which scattering takes place. In the 
case of crystalline silicon (point group Oh), we can define x, y, and z as the 
[100], [010], and [001] directions of the material crystal. The Raman tensor 
for phonon polarization in each of these directions is |22[ [26] 

/O 0 0\ /O 0 d\ /O d 0\ 

R,j,= 0 0 d ,Ry= 0 0 0 ,Ri= d 0 0 , (16) 

\0 d oj \d 0 0/ \0 0 0/ 

where d is the only independent Raman tensor element for crystalline silicon 
since the three optical phonon modes in the T point are triply degenerate 
j22j . The triply degenerate feature of the optical phonons manifests itself 
in the common Raman tensor element d for both LO and TO modes. Note 
that for a polar material the Raman tensor element d for longitudinal optical 
(LO) and transverse optical (TO) phonons will be different due to elasto- 
optic effects. Let us further consider an incident photon with polarization in 
vector direction e^, scattered into polarization ig in this crystal. The Raman 
tensor element in Eqns. [^and[^is given by 

\R\‘^ = \ei-R^-es\\ (17) 

where ^ is the phonon polarization unit vector, and R^ is the corresponding 
linear combination of the Rx, Ry, and Rz matrices [22]. In our calculation 
we fix the incident photon to be traveling in the direction ki = x, having 
polarization di = z with respect to the material crystal. We then quantify 
light scattering for all possible scattered momentum vectors kg. For each kg 
two mutually perpendicular photon polarizations (e^) are considered. Con¬ 
servation of momentum enforces the relationship ki ziz q = kg where q is the 
phonon wave vector. For each q we must consider three phonon polariza¬ 
tions, composed of one ^lo parallel to q and two ^to orthogonal to q . R^ is 
evaluated based on these choices, and the value of the Raman tensor element 
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Figure 2: Band structure and DoS of a diamond structure photonic crystal con¬ 
sisting of air spheres in silicon. The refractive index of silicon is n = 3.44 in the 
calculation, the radius of the air-spheres is r = 0.3a, where a is the lattice constant 
of the photonic crystal. The yellow shaded region denotes the range of the photonic 
band gap. The frequency is in non-dimensional units wa/27rc. 


can be computed. After summing for the two possible scattered pho¬ 
ton polarizations, we obtain the total contributions from LO mode and two 
TO modes to the integral \R\‘^d^l as 6.28|dp and 10.47|dp, respectively. 
These values are valid for all materials in the same crystal structure with 
this particular scattering geometry (x incident light with z polarization). 
For silicon, since d for all three degenerate modes are the same, we collect 
all the contributions in to one single scattering efficiency. This yields a total 
value of 16.76|(ip to the scattering efficiency from all three phonon modes 
in an isotropic DoS medium. 


The net cooling condition (Eqn. 14), demands minimization of the Stokes 
efficiency, thus a photonic DoS design that can fully suppress Stokes scatter¬ 
ing is strongly desired. This is available in various three dimensional pho¬ 
tonic crystal structures, including diamond lattice m, Yablonovite [28], 
woodpile crystal M SO], inverse-opals m. and two dimensional crystal 
stacks [32]. For this study, we select spherical air inclusions in silicon, struc¬ 
tured as a diamond lattice. 

In Fig. [^we present the well-known photonic band structure and DoS of 
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Figure 3: Cooling all Raman-active phonons in silicon: (a) Raman efficien¬ 
cies per unit length Sas and Ss (in intrinsic crystalline silicon patterned with the 
Fig. diamond photonic crystal) are functions of pump frequency. The contribu¬ 
tions from all three Raman-active phonon modes are included. The yellow shaded 
region denotes the band gap for the pump. Band-edge absorption, not included, 
will become significant at pump frequency higher than 260 THz (~ band gap en¬ 
ergy). (b) Calculated cooling and heating efficiency per unit length. The Stokes 
process is suppressed over a wide range resulting in net phonon energy removal and 
the simultaneous cooling of all Raman-active phonons. 


a diamond structure silicon-air photonic crystal, with radius of air spheres 
r = 0.3a, where a is the lattice constant. These calculations are performed 
using MIT Photonic-Bands package m- Since the first Brillouin zone of a 
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diamond structure (face centered cubic) crystal is a truncated octahedron, 
which is nearly spherical [32j . we can make the simplifying assumption that 
the photonic DoS is approximately isotropic. We can then compute the 
photonic DoS through [3l| 


D{u:) = Y, 


FBZ 


dk 

(27r)‘ 



(18) 


where i is the band index, and the integral is over the first Brillouin zone 
(FBZ) of the photonic crystal. 

We now impose this DoS onto the Raman efficiency formulae and assume 
that the escape efficiency for the scattered photons is 1, i.e. all the scattered 
light escapes the medium [35]. In Figure]^ we plot the total scattering 
efficiency per unit length of the Stokes (Ss) and anti-Stokes {Sas) processes 
as functions of pump frequency for the selected photonic crystal. Note that 
we include the contribution of the triply degenerate LO/TO phonon modes 
of silicon into a single Ss and Sas- We invoke the known properties of 
silicon |2T], M = 28.09 amu, N = 2.5 x 10^^ cm“^, \d\ = 1.9 x 10“^® cm^, 
and a = 850 nm as the lattice constant of the diamond structure photonic 
crystal to complete our calculations. Several notes can be made from this 
computation. Pump light lies within the photonic band gap in the 160 - 
203 THz (1476 - 1873 nm) range. Therefore, little to no scattering will 
occur. Pumping in the 203 - 220 THz (1362 ~ 1476 nm) range, however, 
completely suppresses the Stokes process. Simultaneously, the anti-Stokes 
process is enhanced by a factor of ~ 15. 


In Fig. we plot the cooling and heating efficiencies per unit length, 
i.e. ^^^Sas and ^Ss as a function of pump frequency. We observe that 
net cooling can occur within the Stokes suppression range (203-220 THz, 
1362 nm - 1476 nm) as long as Eqn. 14 can be satisfied. The majority of 
optical absorption in crystalline silicon (band gap energy « 1.12 eV) over 
this frequency range comes from three- and four-phonon assisted absorption 
|25j . and is very small. To evaluate the possibility of cooling, we invoke 
experimentally measured absorption coefficients from [231 and compare 
them against the net cooling efficiency per unit length ^j^Sas ~ (pre¬ 
sented in Fig. 1^. As expected, the cooling condition is satisfied for optical 
pumping around 203 to 210 THz (1427 nm - 1476 nm). This demonstrates 
that the total phonon energy removal rate can exceed the energy absorbed, 
leading to net energy removal from the system. 
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Figure 4: Net Raman cooling is achievable in silicon: The net cooling 
efficiency per unit length {-^^Sas ~ can exceed the absorption coefficient 

for the design presented. The absorption coefficient for silicon is extracted from 
[531 [21]. Near 210 THz (1427 nm), the net cooling over comes the absorption. 


5 Practical considerations 

There are a few experimental notes that we must discuss in this context. The 
escape efficiency of the scattered photons, which is the escape probability of 
anti-Stokes photons from the material after scattering, must be engineered 
to approach 100%. Surface states and band edge states are also important 
since they create additional uncontrolled scattering and absorption in the 
material. Several different photonic crystal architectures are known to offer 
the requisite three-dimensional band gaps, including woodpile structure [29l 
SO], inverse-opals m and two-dimensional crystal stacks |32] . The selection 
of a particular crystal design can be made based on our ability to control 
the aformentioned parasitic effects in each structure. 

Cooling efficiency, or heat lift per watt of pump power, is a key metric of 
interest. For photoluminescence based refrigeration |5], the single-pass cool¬ 
ing efficiency is typically around 0.1% per interacting (i.e. fully absorbed) 
photon. In comparison, Raman cooling efficiency is much higher at 5%-10% 
per interacting (i.e. scattered) photon, calculated via the ratio of phonon en- 
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ergy to photon energy. However, since the interaction probability in Raman 
scattering is small, the single-pass net cooling efficiency is very low around 
10 “^ cm“^ (Fig- 0 - To boost efficiency, photoluminescence coolers employ 
pump photon recycling through a multi-pass cavity [5]. However, due to 
the high pump absorption which is required for efficient fluorescence, one 
can only reach a photon recycling factor of about 10, bringing net cooling 
efficiency into the 1-3% regime. In the case of Raman cooling, the cooling 
efficiency could be boosted up by several orders-of-magnitude into the 1% 
regime by adding a high-finesse (Finesse ~ 10^ — 10®) resonator around the 
engineered crystal. This is an advantage of working in the highly transpar¬ 
ent regime of the material. Even further efficiency enhancement can come 
from more careful design of the photonic DoS to express strategically placed 
anti-Stokes resonances. 

6 Conclusions 

In this study, we demonstrate that the spontaneous Raman cooling of trans¬ 
parent semiconductors is achievable when the photonic DoS of the material is 
engineered. Cooling is reached by significantly enhancing anti-Stokes Raman 
scattering over the optical absorption probability, along with simultaneous 
rejection of Stokes scattering events. As a particular example, we show that 
our method can be applied to the net laser cooling of undoped silicon, which 
is currently unattainable by any other laser-cooling method. Furthermore, 
the cooling of all three triply degenerate modes (1 LO and 2 TO) in silicon 
is possible simultaneously using a single telecom wavelength pump. 

In a broader context, Raman scattering and Brillouin scattering are sim¬ 
ilar processes differing primarily by the phonon populations they influence. 
The natural occurence of suitable photonic DoS has already proven to al¬ 
low Brillouin cooling m, thus Raman cooling remains an exciting prospect. 
The ability to exert such control on elementary photon-phonon scattering 
processes through photonic DoS engineering allows us to operate in highly 
transparent regimes far away from any material absorption, band edges, exci- 
tons, or atomic transitions in nearly any semiconductor. Such a broad-based 
method for laser cooling based on Raman scattering can greatly impact our 
ability to control the thermal states of matter. 
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